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All three-loop on-shell HQET propagator integrals n2 n2 

with a loop of a massive quark can be reduced to a set 
of master integrals Q. A reduction algorithm has been 
constructed by solving the integration-by-parts Q re- 
lations using a Grobner bases technique Q . Some mas- 
ter integrals are known exactly, for the remaining ones 
a few terms of their respective e expansion have been 
calculated in Ref. In particular, for the integral 
I2, the 1/e and 0{e^) terms are known. However, in 

some applications more terms of its e expansion are needed (A. Czarnecki, A. Pak, work in 
progress). Here, we shall obtain an exact result for I2 and additional terms of its expansion 
in e. 

We shall consider a more general integral (Fig. |l|): 



Figure 1: The integral /„j. 
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using a method similar to the one used in Q for calculating a three-loop vacuum integral 
Using the HQET Feynman parametrization, we write the one-loop subdiagram as 
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Then, after the Wick rotation pQ = ipEo, we can calculate the integral in d'^^^p: 
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Here Ke I^^^^{ipEo) is an even function oi pEo', Im/^^^^l^Pso) is an odd function, and does 
not contribute to the integral. 

The integral (0) at < can be expressed via the hypergeometric function 2F1: 
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(a similar expression has been derived in Ref. j^). Its square can be expressed via 3F2, due 
to the Clausen identity. In order to perform an analytical continuation to po = ipEO, we 
re-express this 3F2 via three 3F2 functions of the inverse argument: 
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with /t(x) = xr'^{x). This result simphfies at ni = 1: 
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Calculating the integral 
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term-by-term, we obtain 
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In particular, 
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For example, let us consider the convergent integral /i22- Using the reduction procedure 
of Ref. we can relate it to the master integral I2 = Im- 
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From Eq. (11) we have 
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There are several methods to expand the hypergeometric functions in e [§> |^ §]• We 
first follow the method of Ref. . First we express F functions in the hypergeometric series 
as exponents containing 5-sums 
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Some indices are integer at e — s- and some are half-integer; therefore, we get both S'(n; k; 1) 
and S{2n; k; 1). We can re-express 5(n; k; 1) via S{2n; k; ±1) by inserting (1 + (— l)"i)/2 
under the summation sign: 
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Products of 5-sums with upper limit 2n are expressed in terms of single S-sums by means 



of the well-known algebra ||S|, 10, |ll|]. We expand the rational factors in Eq. (|15| ) in e, and 
then expand them into partial fractions. After that, separate sums may diverge, and we 
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introduce an upper limit instead of oo in the outermost sum. Sums with can be 
re-written as sums in n to 2N by inserting (1 + ( — 1)")/2, and those with l/(2n + 1)*^ — 
by inserting (1 — (— 1)")/2: 



^ 1 



n=l 



S{2n;ki, . . . ,ki;xi, . . . ,xi) 



2'=~^[5(2iV; k,h,...,ki; ...,xi) + 5(2iV; k,k^, . . . ,ki; 
:S{2n]ki, . . . ,ki;xi, . . . ,xi) 



^ 1 
^ (2n + l)' 

n=l 



-[5(2iV;A;,A;i, 



. . . - 5(2iV; k,ki, ...,ki; -l,xi, ...,xi) 
- S{2N;k + ki,k2, ■ ■ ■ ,ki;xi,X2, ■ ■ ■ ,xi) + S{2N;k + ki,k2, ... ,ki; -xi,X2, ...,xi) 



(terms vanishing at — > oo are omitted here). After that, all sums divergent at A^ — > oo 
cancel; in the remaining sums, we may set A^ = oo. They are related to the Euler-Zagier 
sums |, |10|. 

Yet another method of expanding certain classes of hypergeometric functions about 
half- integer parameters is provided by the algorithm of Ref. Q , by means of which we have 
calculated the expansion of Eq. (13) up to order 0{e'^). After factoring out an appropriate 
combination of pre- factors, the expansion reads 

^2 



-^122 



TT 



r3(i + e) 



3(1 + 6e) 



1 - vr^e^ + 48,^3e3 



387r^ 
^5" 



- 48 {-K^s - SOCs) £^ + 1152C; 



47937r' 
945 



+ (^39312C7 - ^^^^ - 14407r2C5) + Oie') 



(16) 



In order to check the correctness of this expansion, we have converted the hypergeometric 
functions in Eq. (13) to single Mellin-Barnes representations and subsequently obtained 
the coefficients numerically |l^. We find agreement to at least 14 decimal digits. 

It is now straightforward to obtain the expansion of the master integral I2 from 
Eq. (0): 
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The first two terms have been obtained in Ref. ||T| by a completely different method. The 
next one has been found by A. Pak from the requirement of cancellation of 1/e poles in a 
physical calculation. The last term is then required to calculate the corresponding finite 
part. Yet higher terms in the expansion of I2 can be obtained from Eq. (|T^). 

As a closing remark we would like to mention that the expansion in Eq. ([l6|) is found 
to agree up to order 0{e^) with the expansion of the simple formula 
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We do not know any analytical proof of this result, but high-precision numerical tests 
of Eqs. ( [T3| ) and (|^) for various values of e on the real axis and in the complex plane 
strengthen the conjecture that this result might indeed be valid to all orders. However, 
until an analytical proof of Eq. (18) is found, we are guaranteed to obtain the correct 
expansion only by the analytically derived expression (p!^). 
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